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TO ALL WHOM IT MAY CONCERN: 

Be it known that I, Torsten Franke, a citizen of Germany, whose post 
office address is Reichenbachstr. 34, 80469 Muenchen, Germany, have invented an 
improvement in 

METHOD AND COMPUTER SYSTEM FOR CARRYING OUT A SIMULATION 
of which the following is a 

SPECIFICATION 

FIELD OF THE INVENTION 
[0001] The invention relates generally to a method and system for carrying out a 

simulation and more particularly relates to a system and method of integrating various 
simulation tools. 

BACKGROUND OF THE INVENTION 
[0002] Various computer programs for carrying out simulations are known in the 

prior art. Many of such programs are based on the finite elements method. Such 
computer programs are also referred to as CAD, CAE, EDA and EDM tools. Providers 
of such tools include, for example, the companies ANSYS, INC. (www.ansvs.com) . 
CADENCE DESIGN SYSTEMS, INC. (www.cadence.com ), FLOMERICS, INC. 
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( www.flowmerics .com ), ITI GmbH Dresden (www.iti.de), SimLab Software GmbH 
(www.simlab.de) and The Mathworks, Inc. (www.mathworks.com ). 
[0003] Such computer programs are used for simulation and design purposes in 

all technical areas including for the simulation of technical components and systems. For 
example, the software of the company ITI GmbH allows the modeling of various systems 
such as electrical-mechanical drive trains, internal-combustion and hybrid drives in motor 
vehicles, ships, transporting equipment, machines, gear and brake systems, belt 
conveying installations; the modeling of mechatronic systems, feed drives in machine 
tools, positioning systems, converter lockup systems; the controller optimization of 
hydraulic drive systems, pump and winch drives, valve modeling and activation and also 
the simulation of brake or clutch actuating systems and testbed simulation. 
[0004] One particularly important area for the use of simulation tools is the area 

of electrical engineering and electronics. In particular, such tools are frequently used 
during the design of electrotechnical and electronic equipment and circuits. For the 
simulation of different aspects of the device or circuit concerned several different 
simulation tools are often used. For example, one simulation tool is required to simulate 
the thermal heating caused by the power loss of currents in a device while another 
simulation tool is required for the simulation of its mechanical properties. 
[0005] One disadvantage of the simulation tools known fi*om the prior art is that, 

without being integrated with one another, they address only a partial aspect with regard 
to the simulation and design process. 
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[0006] It has been attempted by what are known as co-simulations to describe the 

problems to be simulated in a simpUfied manner on a common platform (Schwarz, P.: 
KOSIM - A multi-level simulator for the design of integrated circuits. In: H. Reichl (Ed.): 
Proc. MICRO SYSTEM Technologies, Sep. 10-13, Springer, Berlin 1990, pp. 85-90). 
[0007] hi co-simulation, behavioral models are processed. However, the creation 

of such behavioral models requires knowledge of the basic behavior of the object to be 
described (cf Volklein, F.; Zetterer, T.: Einfuhrung in die Mikrosystemtechnik. 
Grundlagen und Praxisbeispiele [Introduction to microsystem technology. Principles and 
practical examples], Vieweg-Verlag, Wiesbaden 2000, pp. 292-326; Schwarz, P.; Haase, 
J.: Behavioural Modeling of Complex Heterogeneous Microsystems; Proc. 1st Intern. 
Forum on Design Languages (FDU 98), Lausanne, Sept. 1998, vol. 2, 53-62; Schwarz, P.: 
Modellierung und Simulation heterogener Systeme [Modeling and simulation of 
heterogeneous systems] set of transparencies; Dresden conference on "Schaltungs und 
Systemsimulation" [switching and system simulation] DASS' 98, Dresden, 05.25.1998; 
Franke, T.; Frohler, U.: Thermal Modelling of Semiconductor Packages; Proc. 
Therminic^OO Workshop, Budapest, September 25-27th, 2000). 

[0008] One disadvantage of known simulation methods with co-simulators, is that 

the automated routines are restricted with respect to their apphcation. As a result, the 
creation of generally usable library elements is not possible on account of the restricted 
areas of vaUdity of these models and a restricted locational resolution. The application of 
these simulation tools, and their linking up, is therefore restricted to selected application. 
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[0009] Further examples of known simulation tools include fluid simulators, 

thermal simulators, switching simulators, simulators for investigating electromagnetic 
compatibility (EMC), mechanical simulators and further CAx, EDA and EDM tools. 

OBJECTS AND SUMMARY OF THE INVENTION 
[0010] It is an object of the present invention to provide an improved method and 

system for simulation. 

[0011] The problem on which the invention is based is solved by the features of 

the independent patent claims. Preferred embodiments are specified in the dependent 
patent claims. 

[0012] The invention allows the results obtained by various simulation tools to be 

stored as library elements in a database. Individual library elements can then be accessed 
according to the question to be simulated, and the library elements can be combined with 
one another to create a simulation model with cross-domain applicability. 
[0013] The system models or parts thereof, obtained in this way can in turn be 

transferred into permanent or temporary library elements. This allows the degree of 
abstraction for definable system areas to be changed and the size of the matrix to undergo 
computation to be reduced without loss in accuracy. 

[0014] In a preferred embodiment, in the operation of bringing the library 

elements together, known as the merger, apart from performing the computation tasks for 
the purposes of the simulation, the result matrix is expanded and graphics are generated 
for the graphic presentation of the simulation result. 
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[0015] It is of particular advantage in this case that a simplified simultaneous 

computation of different physical domains or entities of systems of any desired 
complexity is made possible without loss of resolution or information. 
[0016] Furthermore, the invention provides a link with a development database, 

in which the library elements are stored together with the general properties of the objects 
and which the various simulation computer programs can also access. In this case, the 
elements and objects may be stored in the database directly or as attributive references to 
other sources accessible by computer (for example via a hyperlink). The database can 
also be made available to other tools - such as for example other data ports, simulators 
and devices for machine control. This allows an entirely computer-aided development 
process to be created, with the inclusion of the various stages of the value creation chain 
(vertical integration) and what is known as concurrent engineering (horizontal 
integration). 

[0017] Furthermore, the standardized and computationally reduced model 

structure allows complete integration in the technical systems themselves, such as for 
example in a machine control system. 

[0018] The variable degree of abstraction of the system models generated and 

linked together from library elements is of particular advantage. The linking of the 
various simulation programs via a database system with corresponding library elements 
allows the various simulation programs to be brought together onto a common platform 
while fully preserving the benefits of the individual special simulation programs. This 
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allows incompatible interfaces between the various special simulation programs to be 
eliminated. 

[0019] It is also possible to use different simulation programs of the same 

simulation tool. In this case, it is of particular advantage that the overall amount of 
computation can be reduced by the use of library elements as compared with the 
conventional use of the tool. This is attributable to the library elements or the systems of 
equations reproduced in them being dimensionally reduced in comparison with the 
conventional use of a standard tool. 

[0020] Furthermore, a simulation model created in accordance with the invention 

may itself again be the basis for the application of the method according to the invention 
for creating a further-abstracted system model. 

[0021] A further advantage is the minimization of the overall amount of 

computation which is brought about by the possibility of removing irrelevant matrix 
elements from the matrices containing the transfer functions. This simplifies the 
simulation of complex systems. This advantage exists both in the case of the use of a 
single matrix or library element and when two or more matrices or library elements are 
linked up to form a system model. 

[0022] In particular, the reduction in the amount of computation also allows the 

models to be linked up with automation systems, so that a machine model for the 
purposes of open-loop and closed-loop machine control results directly from the design 
process. 
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[0023] Furthermore, the invention also allows collaboration between remote 

users, such as over the Internet. For example, simulation tasks can be distributed in a 
specifically directed manner along the value creation chain. For instance, the creation of 
a library element may take place at the premises of a supplier of a specific unit. The 
customer can then include this library element in the simulation database, in order to 
generate flexible system models automatically, according to the relevant question, using 
the library elements made available by the supplier. 

BRIEF DESCRIPTION OF THE FIGURES 
[0024] The invention is explained in more detail below on the basis of a preferred 

exemplary embodiment. In the drawing: 

Figure 1 is a flow diagram of an exemplary embodiment of a method 
according to the invention; 

Figure 2 is a flow diagram of an exemplary embodiment for generating a 
system model from library elements; 

Figure 3 is a pictorial diagram showing the linking up of two vector 

models; 

Figure 4 is a flow diagram of a simulation system; and 
Figure 5 is a simplified block diagram showing the incorporation of 
various simulators into the simulation system. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
[0025] Figure 1 is a flow diagram illustrating the creation of various library 

elements for the purposes of simulation. In step 1 , first of all the nodes // and v of a 
vector model are defined with regard to the simulation of a relationship between an entity 
El and an entity E2. 

[0026] The term entity is understood hereafter to mean a specific variable, for 

example a physical, technical or some other variable. For instance, in the case of electric 
current and electric voltage, two different entities are concerned; entities may also 
concern temperature, thermal expansion, mechanical force, mechanical stress, fluidity, 
flow rate, frictional losses etc. 

[0027] The node // concerns those nodes of the vector model which may be 

subjected to an entity El. The node v of the vector model, on the other hand, concerns 
nodes at which a resultant entity E2, for example the vohage, is output. In the case of the 
node V, a distinction is made between two subsets of nodes v* and v"^. 
[0028] The node concerns those nodes v which not only may have the output 

variable of the entity E2 taking place at the node in question but may also be subjected to 
the entity El . On the other hand, the nodes cannot be subjected to the entity El; the 
nodes v"^ serve merely for the output of the output variable of the entity E2 caused by the 
nodes ju being subjected to the entity El . 

[0029] In step 2, the selection takes place of any desired reference point of the 

vector model, which serves as a reference point for the system matrix and for the 
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weighting of the response functions in the later simulation. In step 3, the simulation of 
the transfer functions K^^y then takes place. This involves, for example, selection of tv^o 
nodes v,v and determination of the simulated transfer function with regard to this pair of 
nodes by means of a suitable simulation program. This takes place for all the possible 
pairs of nodes of the vector model defined in step 1 . 

[0030] In step 4, the transfer functions K^^ obtained in step 3 are combined in 

matrices Mej(Kh,,*) and Mei(K^,.+). The matrices Mei(Kh^^.) and Mei(Km^+) thus describe 
the transfer behavior of the system with regard to being subjected to the entity El as the 
input variable and the entity E2 as the output variable. For further computations, the 
matrices may also be represented as scatter matrices, chain matrices or in some other 
form. Corresponding representations are known from Gerstenmaier et al. Proceedings of 
the 4th International Conference on Modeling and Simulation of Microsystems 
(MSM2001), March 19-21, 2001, pp. 84-87. The first-mentioned matrix in this case 
comprises those transfer functions Km which result from pairs of nodes of nodes \i and 
nodes v'*'. On the other hand, the second matrix Mei(Km^+) comprises those transfer 
functions K^^ which result from pairings of nodes \i and nodes v"". The two matrices 
generated in this way are then stored in a database. 

[0031] In step 5, the definition of a further vector model with regard to the effect 

of the entity E2 on an entity E3 then takes place. The definition of the vector model takes 
place in step 5 in a way analogous to the procedure in step 1 with regard to the nodes ju^ 
and v'. In a corresponding way, the selection of a further desired reference point takes 
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place in step 6 and the simulation of the corresponding transfer functions K^'y takes 
place in step 7. For the simulation of the transfer functions in step 7, a further specialized 
simulation program, which can simulate the relationship between the entity E2 and the 
entity E3, is used. 

[0032] In a way analogous to step 4, the output of the corresponding matrices 

Me2(Kjli \v^') and Me2(Kju \ V+ ') then takes place in step 8. 

[0033] Figure 2 illustrates the steps for combining the library elements of the 

database generated in steps 4 and 8 of Figure 1. In step 20, the linking up of the vector 
models takes place. This is performed by assignment of the corresponding input and 
output variables. For example, a subset of the nodes v may serve as input variables for 
the nodes //' of the second vector model with regard to the mapping of the entity E2 onto 
the entity E3 (cf Figure 1, steps 1 and 5). 

[0034] In a way corresponding to the linking up of the vector models in step 20, a 

linking up of the corresponding matrices takes place in step 21. This linking up of the 
matrices results in step 22 in a system model for the simulation of the relationship 
between the entity El and the entity E3. In this case, it is particularly advantageous that 
this system model can be generated automatically from the previously generated library 
elements, without further simulations being required. 

[0035] Figure 3 shows a vector model 30 with a reference point and nodes //, 

v"^ and v"^. The arrows in the connections between the reference point and the 
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corresponding nodes define the direction from which they are subjected to the entity El. 
By definition, the nodes v"^ are not subjected to the entity El. 

[0036] Furthermore, Figure 3 shows the Unking up of the vector model 30 with a 

further vector model 31. The vector model 30 in this case corresponds for example to the 
vector model which was defined in step 1 of Figure 1. The vector model 31 corresponds 
for example to the vector model defined in step 5 of Figure L In the example considered, 
the linkup of the vector models takes place via their reference points and P2. 
[0037] The generation of the matrices or library elements, for example, the 

matrices ME](KM^y^) and MEl(Kpi^v-\-) in steps 4 and 8 of Figure 1, may in this case take 
place by a node // being subjected to a unit step Freference-i^^ The transfer function Kfi^y 
can then be determined from the response AD to the unit step, determined between the 
point V and the reference point Pj or P2. 

K (t) = AD, U)/FZ (0 I^A^ 
(/<0) = 0 ; F^, > 0) = const ;t 0 

[0038] For each pair of nodes, such a transfer function K^^y is separately 

determined, so that altogether a number of // v matrix elements is obtained. 
[0039] One possible computation routine for calculating the output variables and 

the nodes by means of boundary conditions Ky^^equilibrium(0 the transfer 
functions when induced by a function Fy^(t) is given below: 
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(0 = K,...^„,,^^ (/) (/) + (0 I V V * 



[0040] With regard to the formulas given above, it is presupposed that the 

reference point or P2 does not lie at a node which is assigned to a boundary condition 
^v'^-equilibriumO) ' principle, however, an entirely free choice of the reference point is 
possible. With a free choice of the reference point and with the nodes // subjected to 
feeding, the following extended system of equations is obtained: 



(changed indexing for F(t)) 
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Vi/* (in addition) 



(in addition) 



with Kreference^egudibnumft) " boiuidary conditioii at the reference node 

^reference^^^ ; F^^(t) - component of the inducing function Fre/erenceft) or Fy^(t) which 
is dependent on the respective boundary condition Kreference-equihbnum(t) or Kv*,equihbnum(0 
Preference ' Component of the inducing function Freference(t) or Fy^(t) which is 

independent of the respective boundary condition Kreference-eguihbnum(t) or Ky*-equiiibnum(t)^ 

[0041] For this calculation, only * v* = instead of * v matrix elements are 

required. The remaining |a(v-v*) elements do not contribute to the actual solution and 
can therefore be removed from the matrix. 

[0042] An example of an expansion routine for the determination of the elements 

not to be considered in the calculation is given below: 
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[0043] Furthermore, an expansion routine for the determination of variables not 

explicitly considered in the model, with reference back to known object properties, is 
given below: 

[0044] In this case, the gradient is a direct result of the above calculation and is 

an object attribute which can be called up from the database. 

[0045] Figure 4 shows a flow diagram of an integrated system for the creation of 

library elements, for the creation of system models, for carrying out simulations and for 
the output of simulation results, which integrates various simulation programs with one 
another. 

[0046] In step 40, the "merger", that is to say the integrated simulation method, is 

called up. In step 41, it is determined whether the program is to be closed again. If this is 
the case, the closing of the program takes place in step 42. 

[0047] In steps 43, 44, 45 and 46, the desired operating mode of the program is 

then called up. In step 43, the "create element" mode can be chosen for the creation of a 
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new library element. In step 44, the "create system" mode, in which a system model is 
created and stored in the library, can be chosen. 

[0048] In step 45, the "calculate system" mode, in which a simulation using a 

system model is loaded from the library and the output of the result being stored into the 
library takes place, can be selected. In step 46, the "present results" mode can be chosen. 
In this mode, the results are loaded from the library and output in the desired format. 
[0049] If the "create element" mode is chosen in step 43, the corresponding tool is 

called up in step 47 and the corresponding vector model is loaded in step 48. In step 49, 
the transfer functions for the corresponding model are determined by way of a 
corresponding specialized simulation program and are extracted as parameters. 
[0050] In step 50, the tool is closed again and, in step 51, the library element is 

initialized. The library element is then stored in a database library in step 52. 
[0051] If the "create system" mode is chosen in step 44, the corresponding 

database library is opened in step 53 and the required library elements or matrices are 
loaded in step 54. In step 55, these elements are arranged for linking up with one other, 
so that an initialization of the system model can take place in step 56. This system model 
is stored in the database library in step 57. 

[0052] If the "calculate system" mode is chosen in step 54, the corresponding 

database library is again opened in step 58, in order to load the previously generated 
system model in step 59. 
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[0053] In step 60, the calculation for the simulation by means of the system 

model is initialized and the corresponding calculation is performed in step 61. The 
results of the calculation are stored in the database library in step 62. 
[0054] If the "present resiilts" mode is chosen in step 46, the database library is 

again opened in step 63, in order again to load the results previously stored in step 62. 
Furthermore, the object properties are loaded in step 65, in order to carry out an extension 
of the results in step 66. In step 67, the user can select the type of presentation of the 
output or the presentation format. In step 68, the results are then output in the desired 
presentation. 

[0055] To simplify use of the integrated simulation system further, it can also be 

made accessible via a browser over the Internet. This arrangement facilitates the 
implementation of linked-up applications. 

[0056] If a database library with Ubrary elements generated in advance already 

exists, individual functionalities (create element, create system, calculate system, present 
results) can also be used independently and/or in parallel. 

[0057] Figure 5 shows a block diagram of the connection of the merger tool of 

Figure 4 to various simulation programs. The system of Figure 5 has a general 
development database 70, in which various data for the purposes of the development of a 
technical system or a technical component are stored. Textual information, drawings and 
further related data can also be stored in the database 70. 

[0058] The simulation programs 71 to 76 including the fluid simulator 71, 

thermal simulator 72, switching simulator 73, EMC simulator 74, mechatronic simulator 
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75 and further CAx, EDA and EDM tools 76, access the general development 
database 70. These simulation programs 71 to 76 are likewise connected to the merger 
tool 77. 

[0059] This arrangement creates an integrated development environment, which 

allows even unskilled users to perform complex simulation tasks. The output of the 
merger tool 77 may in this case be used to load a machine model generated by the merger 
tool 77 into the machine control system 78. Furthermore, the merger tool 77 is also 
connected to the general development database 70, in order to store the results of the 
merger tool 77 in the general development database 70. 
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